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We investigate two-dimensional turbulence within the Instanton formalism which determines the 
most probable field in a stochastic classical field theory starting from the Martin-Siggia-Rose path 
integral. We perform an approximate analysis of these equations based on a variational ansatz using 
elliptical vortices. The result are evolution equations for the positions and the shapes of the vortices. 
We solve these ordinary differential equations numerically. The extremal action for the two-point 
statistics is determined by the merging of two elliptical vortices. We discuss the relationship of this 
dynamical system to the inverse cascade process of two-dimensional turbulence. 
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I. INTRODUCTION 

Over one and a half centuries after the discovery of 
the Navier-Stokes equations, the theoretical understand- 
ing of turbulent flows remains quite limited although phe- 
nomenological theories (see Jj) capture the main charac- 
teristics of turbulent flows quite well. However, no gener- 
ally accepted derivation of the phenomenological theories 
from the basic hydrodynamical equations has been possi- 
ble. The challenge consists in connecting dynamical sys- 
tems theory with nonequilibrium statistical physics [2j. 

The study of two-dimensional turbulence is motivated 
by the fact that whenever one spatial direction is signif- 
icantly constrained (e.g. layers in the atmosphere or all 
kinds of surface dynamics) the flow can be described as 
quasi two-dimensional. In a two-dimensional space, the 
behavior of turbulent flows differs significantly from the 
three-dimensional case. These differences arise due to an 
additional conservation law for the enstrophy which leads 
to a new phenomenon called the inverse cascade discov- 
ered in 1967 by Kraichnan [3]. Numerically, a stationary 
inverse energy cascade can be generated by a stochastic 
forcing of the flow at small scales while dissipating energy 
at large scales. 

The stochastically forced vorticity equation of two- 
dimensional turbulence can be treated as a stochastic 
classical field theory using the Martin-Siggia-Rose path 
integral formalism ^ . Path integrals may be treated ei- 
ther perturbatively using renormalized perturbation the- 
ory or within the Instanton formalism by an expansion 
of the field around the most probable field configuration. 
This formalism has been applied by several groups to the 
Burgers equation (see [SHE]), to the problem of turbulent 
advection and to estimate the tails of the velocity prob- 
ability distribution of fully developed turbulence [9] . We 
also refer the readers to earlier works by Falkovich [TD] 
and Moriconi [llj . 
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The purpose of the present contribution is to in- 
vestigate the dynamics of the Instanton equations for 
the inverse cascade in two-dimensional turbulent flows. 
Thereby, we do not focus on the concrete evaluation of 
the extremal action and the shape of the probability dis- 
tribution. As has been pointed out by [3], a reasonable 
estimation of the probability distribution has to take into 
account fluctuations around the Instanton. We rather ad- 
dress the question whether one can identify signatures of 
the mechanism underlying the inverse cascade already on 
the basis of the A^-point Instantons. 

Our approach consists of an ansatz for the Instanton 
fleld based on a superposition of elliptical, deformable 
vortices. We derive the most probable evolution of these 
vortices by extremalizing the Martin-Siggia-Rose ^1 ^tc- 
tion (see eq. This ansatz is motivated by the obser- 

vation that the Instanton equations in the inviscid limit 
and in the presence of small scale forcing are related to 
the dynamics of point vortices (see [T^HIl])- 

We design a numerical procedure to solve the resulting 
evolution equations for the positions and the elliptical 
shapes explicitly for two vortices. This gives the two- 
point statistics in the Instanton approximation. Our aim 
here is not the complete determination of this action, 
which is a project for future research. The aim of the 
present paper is to investigate the deterministic dynamics 
of the Instanton equations. Our variational ansatz using 
elliptical point vortices demonstrates that the action is 
essentially determined by the interaction of two elliptical 
vortices. Thereby, the elongation of the main axis of the 
elliptical vortices plays a crucial role. This elongation 
leads in the case of the interaction of vortices with like- 
signed circulation to an accelerated relative motion (see 
also [H]). 

The structure of this paper is as follows. We start 
with the investigation of the dynamics of Instanton point 
vortices which do not exhibit an inverse cascade. We 
present a generalized model that takes into account the 
deformation of vortices due to strain and calculate the 
interaction of two vortices numerically. We discuss the 
relationship of this model to the inverse cascade process. 
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II. INSTANTON EQUATIONS 

We start from the two-dimensional vorticity equation 
for incompressible fluids with lj = V x m 

dtuj{x, t) + u{x, t) ■ Vcj(a;, t) - z^Aw(a;, t) = /(a;, t) (1) 

where f{x,t) denotes an external Gaussian white noise 
forcing which is specified by 

{f{x,t)f{x',t'))=Qix-x')S{t-t'). (2) 

The velocity field can be expressed in terms of the vor- 
ticity field via the Biot-Savart law according to 



u{x, t) 



(Px' x~ x' 



27r la; — x\ 



uj{x',t). (3) 



Within the Martin-Siggia-Rose formalism [4], all mo- 
ments are derived from the partition functional which 
is given by the path integral 



VujVu} e'^ 

with the extended Martin-Siggia-Rose action 



(4) 



S = I (PxdtLu{x,t) 
i'Auj{x, t) + ■ 



dtOj{x, t) + u{x, t) ■ \7uj{x, t) 



d^x' Q{x~ x')6j{x' ,t) 



i / dixAt'q{x,t)uj{x,t) . (5) 



The partition functional Z can be related to the 
Laplace transform of the A^-point vorticity probabil- 
ity distribution by evaluating Z for the field ri{x,t) = 

The Instanton equations are derived by performing a 
variation of the action ([s]) with respect to the vorticity 
{SS/6uj = 0) and the auxifiary field {6S/66j = 0). This 
procedure yields two coupled partial differential equa- 
tions that describe the evolution of the most probable 
field configuration 

dtuj — — u ■ Vui + vAui — i J d^x' Q{x — x')uj{x' , t) 

dtCb = - 

— u ■ Vuj — vAuj + rj{x, t) . 
The extremal action is then given by 



27r 



x' ^, . e^-K{x-x') (6) 



iS'c: 



d^x' Lij{x,t)Q{x- x')Cj{x' ,t) . (7) 



The equation for ui has to be evaluated forward in time 
from t* < Q io 0, whereas the equation for Cj has to be 



solved backward in time from to t* < 0. This is in 
accordance with the different sign of viscosity. For the 
auxiliary field we have the initial condition 



(a;,0) = ^^a.iS{x - Xi) 



(8) 



4=1 



which results from considering the n-point statistics of 
the field at spatial points Xi at time t = by setting 



III. ZERO- VISCOSITY LIMIT 

In order to motivate our further treatment, let us con- 
sider the evolution equation for the auxiliary field in the 
limit of vanishing viscosity under the neglection of the 
third term on the left-hand side of equation ([T]) 



dtilj + u ■ Vtj = ai5{x — Xi)S(t) 



(9) 



The solution to this equation is given by a sum over delta- 
functions 



Ld{x, t) = ai5{x — X{xi, t)) 



(10) 



i=l 



where X{xi,t) denotes the Lagrangian path (for t < 0) 
of a passive particle moving in the velocity field u{X,t) 
ending at Xi at time t = governed by the evolution 
equation 



^^X{x,t)^uiX{x,t),t) 



(11) 



with the velocity field related to the field uj{x, t) by Biot- 
Savart 's law. 

The temporal evolution for the Instanton field Lo{x,t) 
in the zero- viscosity limit is then given by 



dtio + u ■ Vw = — Q{x — X{xi, t))c 



(12) 



This equation can also be solved using the Lagrangian 
path X(x,t) by evaluating 

ujiX{x,t)) = 

-VaW dt'Q{X{x,t')- X{x„t')). (13) 

For small-scale stirring we can approximate the correla- 
tion Q{r) by a 5-function. 

Finally, we have to determine the velocity field, which 
is obtained by a combination of Biot-Savart 's law and the 
representation of the vorticity field 



u{x, ^) = X! J ^ 27r 



\x — X(xi, t' 



(14) 
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For the sake of simplicity, we approximated Q{x — x') 
by a delta function, q5{x — x'). We obtain an extension 
of the usual point vortex dynamics by the fact that the 
circulations is time dependent oc (i — t*)aiq. 

This dependence of the initial time t* can be taken into 
account by a time shift t = t ~ t* . The Instanton point 
vortex dynamics is given by (xilr) = X{xi, r + <*)) 



X 



(15) 



The transformation r' = ^r^ yields the usual point vor- 
tex dynamics. Therefore, we conclude that the inviscid 
Instanton point vortex dynamics is similar to the usual 
point vortex dynamics in an accelerated timeframe which 
consists of replacing the r-dependence by a depen- 
dence. 

On the basis of the point vortex approximation, there 
is no mechanism which could be related to the inverse 
cascade. 



IV. DEFORMABLE VORTICES 

In this section we sketch a generalization taking into 
account the finite extension of the vortices as well as a 
deformation of the vortices due to shear generated by dis- 
tant vortices. This will motivate our approximate treat- 
ment of the Instanton equations. 

We start with the auxiliary field, which for i — >■ actu- 
ally has to tend to a superposition of (5-functions, i.e. to 
a superposition of point vortices. For negative times, vis- 
cosity will broaden the vortices whereas presence of the 
advective term and the type of vortex stretching term 
will lead to a deformation of the circular shape. Thus, 
the straightforward extension is an ansatz in terms of 
elliptical point vortices 



j{x,t)^j2- 



,-i(x-%(t))cj l'(x-%(t)) 



j=i 27rA/detC, 



(16) 



From the condition ([s]) we conclude that 

fj(0) = aj,Cj{0) = 0, and %(0) = Xj . (17) 

A similar ansatz is performed for the Instanton field 
uj{x, t) 



cj{x, t) 



E 



r,(t) 



-i(x-a;j(t))Cj(t)<"^'(x-Xj(t)) 



^ 27r^detQ(t) 



(18) 



The initial condition for t* — > — oo, u}{x,t) = is defined 
by Tj{t*) — 0, Cj{t*) and the initial locations Xj{t*). 
The initial condition uj{x,t*) = can be realized by 

r,(r) = 0. 



It is evident that we could extend our ansatz by in- 
cluding, besides the elliptical ansatz, higher-order terms. 
Such an expansion with an additional averaging over 
the higher-order contributions is needed in order to go 
beyond the bare Instanton approximation. This pro- 
gram would lead to an effective action in much the 
same way as has been demonstrated by Falkovich et 
al. [3] in their treatment of the statistics of the vor- 
ticity increment of the direct cascade. They separated 
the fields into an axisymmetric and a nonaxisymmetric 
part, where the nonaxisymmetric fluctuations where in- 
tegrated out to obtain an effective action for the axisym- 
metric part, which, in turn was treated in the Instan- 
ton approximation. We conjecture that such an effective 
action could be characterized by a modified noise corre- 
lation Q{x — x') as well as a modified Biot-Savart law, 
u{x,t) = Bz VF{—A)u}{x,t), where F{—A) is a func- 
tion of the Laplacian. 

In the present article we will be content with the el- 
liptic vortex approximation. Our main point here is that 
the behavior of the elliptical vortices depends on the signs 
of Tj and, hence, on the sign of aj. We will explic- 
itly demonstrate that two like-signed vortices will attract 
each other due to this deformation. 



V. ACTION FOR ELLIPTIC VORTICES 

In this section we consider the Martin-Siggia-Rose ac- 
tion for a field which consists of a superposition of N 
elliptical vortices. Each vortex is determined by its lo- 
cation Xi{t), its amplitude Ti{t) and a matrix Ci{t) de- 
termining the orientation and the semi-axis of the el- 
liptical vortices. Similar quantities are introduced for 
the auxiliary field. By introducing these collective vari- 
ables, the Instanton equations, which are partial differen- 
tial equations, are approximated by ordinary differential 
equations, which are the Euler-Lagrange equations for 
the MSR-action, written in collective variables. These 
equations turn out to be extensions of the inviscid point 
vortex Instantons. The inclusion of the elliptical shape 
of the vortices, however, will lead to a modified vortex 
dynamics. 

As a starting point we consider the Martin-Siggia-Rose 
action (eq. (Is])) in Fourier space 



ld?k dtu{-k, i) |5fw(fe, t) + uk^Luik, t) 
ik ■ J d^k'u{k')uj{k - k',t)uj{k\t) 



A^kAtCj{--k,t)Q{k)Cj{k,t) (19) 



where we defined it(fe) = i^^}^. We make an elliptic 
ansatz for the vorticity and the auxiliary field according 
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to 



1=1 



(20) 



Ci and Ci are 2x2 symmetric matrices and Ti and 
are the circulations of each vortex or eUiptic structure in 
the auxiUary field. 

Here, we assume that the quantities Ti{t), Xi{t), Ci(t) 
as well as the corresponding quantities for the auxiliary 
field depend on time. Extremalization of the MSR-action 
will lead us to an approximation of the Instanton equa- 
tions in terms of a coupled set of ordinary differential 
equations. As we have discussed in the previous section, 
we include the deformation of the vortices since this de- 
formation seems to play a major role for the energy trans- 
port in the inverse cascade. 

We have to specify the boundary conditions for the 
fields with respect to time. The conditions for the auxil- 
iary field are easily formulated. From the representation 
dsl) we obtain 



i,(0) = X, and (7,(0) = 0. 



(21) 



The initial condition for the field uj{x,t) has to be for- 
mulated for t — >■ —oo, in particular we have ui{x,t — >■ 
— c») — 0. As we will see below, this can be achieved by 
putting Ci{—oo) = 0, i.e. a infinitely extended vortex. 
Another possibility would be to set Ti{t*) = for a large 
negative time t* . The reason will become more evident 
below, when we formulate the Euler-Lagrange equation 
corresponding to the MSR-action for the elliptical vor- 
tices. 



We plug the ansatz into the action ( 19 1 and integrate 



the action with respect to k which yields 



^ — 1 A — 1 



i=l J=l ■ 

- i^FfeAjfc(VO - vVaV, + ]^^Kj{V,) \w,j , 

k=l ^ J 

(22) 

where denotes Vj,; and we defined 



(23) 



Kj{\/,) = iQ(iV,)e-(*^-"^)-^*+5^'(^^-^^)^' 



W,, = i 



2tt 



'dct[Ci+Cj] 



Equation ( 22 1 is the action for elliptic vortices and elliptic 



structures in the auxiliary field according to the ansatz 
from eq. (20). We will derive the evolution equations 



for the parameters in the ansatz by applying a variation 
according to the Euler-Lagrange equations. To this end, 
we make the following assumptions. 



1. We assume that we have the same number n ~ m ~ 
N of elliptic vortices and elliptic structures in the 
auxiliary field. 



2. The elliptic structures in the auxiliary field are po- 
sitioned such that the i-th vortex is always close to 
the i-th elliptic structure, which means 



\xi{t) — Xi{t) \ ^ Size of vortices . 



(24) 



3. The vortices and elliptic structures are isolated, 
hence their overlap can be omitted. This yields 
Vi ^ j 

\xi{t) — Xj{t)\ ::§> Size of vortices (25a) 
and Vi 7^ j 

\xi{t) — Xj{t) \ ^ Size of elliptic structures . (25b) 



Hence, we only account for the diagonal terms in eq. 



(22). Based on this approximation we pursue the inves- 



tigation of the action 



i=l 

N 



Fj - 1 . . . 

Li 2 



i^ FjA,j(V) - i^VIV + 2 T^^i(V) \Wu , (26) 



where we have set i = j and then renamed k to j. 



VI. EVOLUTION EQUATIONS FOR 
COLLECTIVE COORDINATES 



We formulate the evolution equations for our Instanton 
model of elliptical vortices and the corresponding auxil- 



5 



iary field 



until first order in Ci + C, and define 



N 



Xi ^^TjUi.j{xi - Xj) + q-^{xi - Xi{t)) 



N 



r,, = 



0, 

N 



{q 



Ci — Ci 



(27) 



Xi — ^ ^ C/ij (aji a; j ) , 



AT 

Cj — ^ ^ l^^ij {xi — Xj) Ci + CiSfj {xi — Xj) 



2^1 



in collective coordinates as the main result of our ap- 
proach. Note that we accounted for the previously men- 
tioned assumptions in the derivation of these equations 
which can be found in the appendix [K\ 

These equations allow for a subset of solutions for 
which the locations of the vortices of both fields Xi{t) 
and Xi{t) coincide. We will discuss this case later. 



V{x) 



1 



27rr2 V ^ 



-y 



(29) 



We can write Uij{r) with C^l^'' = i(C^l^ + C^^J) as 



f/.(r-) = V(r)+E^i?^a^^W (30) 

a,/3 — 1 ^ 

where cj^} denotes the component of the matrix Ci. 
For the deformation we have to evaluate 



1+ (C.+Q)V. 



^ V^^r (31) 



which we can write as 



with the matrix 



M{x) = 



1 



2xy 2y^ - r 
'2xy 



(32) 



(33) 



An analogue calculation for the auxiliary field yields 
U,ir)-Vir^^tcti!^l^ir) (34) 

Q,/3 — 1 



VII. TWO-POINT INSTANTON 

We solve the evolution equations explicitly for two vor- 
tices with same-signed circulations. In accordance with 
the assumptions made earlier, we allow for two ellipti- 
cally shaped structures in the auxiliary field which have 
to tend to infinitely small objects once time t tends to 
0. This means that at time the auxiliary field consists 
of a superposition of two delta functions (see eq. ([s])) 
and hence is connected to the two-point statistics of the 
vorticity field. We call the corresponding Instanton the 
two-point Instanton which provides an insight to the in- 
teraction of two vortices. 

We want to derive the dynamics of the two-point In- 
stanton to understand the underlying vortex dynamics. 
To this end, we expand the advection and deformation 
and make use of an iterative calculation scheme to obtain 
the solution for the vorticity and auxiliary field. 

Let us consider the calculations in some more detail. 
We want to calculate the evolution of two elliptic vor- 
tices and two elliptic structures governed by the evolution 
equations (27). 

To this end, we expand the advection for i ^ j |19) 



and 



U^j (r) 



1 + ;^V^ (a-fQ)V. 



2Ti\r\ 



(28) 



S.,,{r) ^ M{r) + 



2 

E 



c. 



{ij) 

Q/3 



_9 d_ 

dxa dxp 



M{r) 



(35) 



~ (ii\ ^ (i\ (j\ 

where C/jj denotes C/h; -I- Cjjil and r Xi ~ Xj. 

The auxiliary field satisfies the initial condition 
u;{x, 0) = J2i=i ctiHx — Xi) according to eq. Q which 
corresponds to two point vortices. The vorticity field is 
fixed at a time t* < 0. Because the initial conditions 
for the two fields have to be fixed at different times, a 
straight forward calculation is not feasible. 

We solve the equations for two vortices and two el- 
liptic structures in the auxiliary field using an iterative 
technique. We start from the well known solution for 
two point vortices (this emerges in the limit of vanish- 
ing Ci and Ci) located at a;i and X2 or Xi and X2 re- 
spectively with Xi{t) = Xi{t) for i = 1,2. The initial 
condition for the auxiliary field consists of two point vor- 
tices located at Xi and X2 so that one can argue that uj 
will behave similar to the motion of a two point vortex 
system. Now we calculate the evolution of the vorticity 
field (jj according to the set of equations ( |A9[ ) in the limit 
Ci — ?> 0. Hence, we have Xi{t) — U{'dt)xi{0) where U 
denotes the two-dimensional rotation matrix 50(2) and 
I? = (Fi -I- T2)/ {2Tr(P). We use the approximations for 
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the advection Uij and deformation Sij from equations 
(30) and (32) to calculate Xi{t) and Ci{t). The ordi- 



Distance 



nary differential equations can be solved with an explicit 
Runge-Kutta method. 

Now we can calculate the auxiliary field w for Ci 
if we plug in the previously calculated solution for the 
vorticity field given by r,;(i), x{t) and Ci{t). Because 
of dtTi « the circulation Ti(t) grows linearly in time 
according to qVi [t — t*) + Ti {t* ) . With the obtained so- 
lution for the elliptic structures in the auxiliary field we 
can calculate a more precise solution of the vorticity field 
and successively repeat these steps. Due to the smallness 
of the coupling parameter q compared to the the circu- 
lations {q ~ 0.05 where we chose the circulations in the 
order of unity) this iterative procedure converges quickly. 

The result for same-signed vortices is shown in figure [T] 
Initially circular shaped vortices are deformed with time, 
forming long elliptic structures. These structures then 
approach each other which is shown in figure [2j where we 
plotted the distance between the centers of the vortices. 



VIII. REDUCED EQUATIONS 



We want to discuss the relationship of the elliptic 



model described by the set of evolution equations (27) 
to the inverse cascade in more detail. To this end, we 
formulate a reduced version of the evolution equations. 
We perform the following approximations: 

1. The position of the center of each vortex is equal 
to the center of each structure in the auxiliary field 





Vorticity field 




Auxiliary field 



Figure 2: Evolution of the distance between the centers of the 
vortices (solid line) and the centers of the elliptic structures in 
the auxiliary field (dashed line). One sees the tendency of the 
vortices to approach each other rapidly once the elongation 
becomes dominant (vertical dotted line). 



at each time t between t* and 0. This means 

x, = x, Vte[r,0] (36) 



and hence the term qY:{xi~Xi{t)) in (27) vanishes. 



2. The structures in the auxiliary field are infinitely 
small and can be described by 5 functions. In terms 
of the collective coordinates, this means 

a = V<e[r,0]. (37) 

As a consequence, the coupling term in the evolu- 
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tion equation for Ci reduces to 



(38) 



and hence the evolution of Gi becomes decoupled. 

We obtain a reduced set of evolution equations for the 
vorticity field that implicitly contains the evolution of 
the auxiliary field given the above approximations. One 
can think of it as an enslavement process where the evolu- 
tion of the point vortex like auxiliary field follows instan- 
taneously the centers of the elliptic vortices. It shall also 
be noted that the approximation Ci{t) = is consistent 
with the initial condition at time t ~ Q. The resulting 
equations for the vorticity field are 

Ti = qti , 

N 



N 

Ci — ^ ^ V j \^^Sij{xi Xj)Ci + CiS^jlxi X 
f, 



(39) 



01 



2vl. 



Given the solution to ( 39 ) the evolution of the auxiliary 
field is well known. The circulation is constant, the 
position Xi is equal to Xi and the matrix Ci is zero. Nev- 
ertheless, the effect of the auxiliary field is still present 

and taken into account by the term q^{Q — Ci). 



It is interesting to mention that ( 39 1 is equivalent to 



the equations for elliptic vortices presented by Friedrich 
and Friedrich in |16j from which they derived the ro- 
tor model. The term q^riQ ~ C'i) can be considered an 
overdamped spring that forces the shape of the elliptic 
vortices to a certain shape. In the limit of strong elonga- 
tion the elliptic vortices can be approximated by pairs of 
point vortices connected via an overdamped spring with 
the equilibrium distance Dq and an appropriate spring 
constant 7/2. This leads to the rotor model which ex- 
hibits an inverse cascade by cluster formation of same- 
signed vortex pairs starting from arbitrary initial condi- 
tions. See [TB] for a detailed discussion. 

In summary, the effect of the auxiliary field yields the 
extension of the point vortex model that is responsible for 
breaking the F^ — > — Fi,i — > —t symmetry giving rise to 
an inverse cascade as demonstrated by the rotor model. 



IX. DISCUSSION 

We investigated the dynamics of the Instanton equa- 
tions for the inverse cascade in two-dimensional turbulent 
flows. Because the dynamics of point vortices as consid- 
ered in the limit of vanishing viscosity do not lead to an 



inverse cascade we decided to introduce a more sophisti- 
cated model taking into account a finite spatial extension 
of the vortices as well as their ability to perform deforma- 
tions on the basis of shear generated by distant vortices. 

To this end, we derived the most probable evolution of 
elliptic vortices by extremalizing the Martin-Siggia-Rose 
action under certain assumptions. The main result is 
given by the evolution equations in collective coordinates 
(p7|. 



We focused on the interaction of two same-signed el- 
liptic vortices where we observed the elongation of the 
main axis which led to an acceleration of the relative 
motion. In our opinion, the combination of the elonga- 
tion and the acceleration of the relative motion which 
leads to a reduction of the distance between two same- 
signed vortices is the signature of the inverse cascade. A 
reduced set of equations for the evolution of the vorticity 
field was derived which leads directly to the rotor model 
|16) . It exhibits an inverse cascade by cluster formation 
of same-signed rotors which correspond to infinitely elon- 
gated elliptic vortices. 

In summary, we investigated the two-point Instanton 
by evaluating the interaction of two deformable same- 
signed vortices. We identified the underlying dynamics 
in accordance with [17] and [18] as elongation and thin- 
ning of elliptical shaped vortices due to the shear and 
hence the acceleration of the relative motion which leads 
to enhanced transport of energy and enstrophy. We con- 
clude, that the Instanton calculation supports the view 
that elongation and thinning of small scale vortices and 
the resulting clustering dynamics lead to the formation 
of the inverse cascade in two-dimensional turbulence. 

It remains an interesting task for the future to gener- 
alize our approach by modifying the ansatz. One could 
think of vortices with deformable main axes so that one 
would expect to obtain S-shaped structures or an exten- 
sion to the exponential of the form u oc exp[(a;C~^a;)"]. 
The evaluation of the two-point statistics within the ellip- 
tical model poses an interesting challenge which consists 
of the evaluation of the transition probability which we 



sketched in the appendix (see eq. (B2)). To this end, one 
has to consider solutions to the Instanton equations that 
deviate from the point vortex vicinity and thus disobey 
the assumptions we made for our approach. 



Appendix A: Evolution equations for elliptic vortices 

We start from the Lagrangian 



N 



F, 



1 



i=l 

N 



i^F,^,,(V) - ^.VIV + ^^i^.(V) \W., , (Al) 



which one obtains immediately from the action ( 26 1 



The action and thus the Lagrangian do not depend on 



8 



the time derivative of the auxihary field so the Euler- 
Lagrange equations for the variation with respect to lo 
reduce to 



dC 



(A2) 



with fi E [f i, Aij. , Ci^J for k,l — (1,2). The variation 
with respect to Ti yields 



N 



= < ^ - • V + ^vav - i^r,A,,(v) 



i^yiV + i^ifz(V) \ Wu , (A3a) 



the variation with respect to Xi leads to 
= -V + ivaV 

N ^ 

- i TjA.j (V) - i/VXV \ V Wu (A3b) 
j=i J 

and the variation with respect to Ci provides 

fr 1 ^ 

- <^ - ±, • V + -vav - i^r,A,,(v) 



i/VIV + i-^K.iV) WV^ VF,, . (A3c) 



We make a Gaussian ansatz for the correlation function 
Q according to Q{k) = qexp[—^kQk]. The width of Q 

determines the forcing scale A. We expand Aij{W) and 



Kiiy) in equations (A3a)-(A3cl according to 



1 + -VC.fc' + -fc'C.V 
2 2 



tt(fc') • V (A4a) 



and 



K,{V) « ig 



1 - - a;,) • V + -V(g + Q - Q)V 



and define {af^'' and is of A:-th order in V) 

N 



a. 



(1) 



J2^jU,j{x, - ajj) • V, 



(2) 



where Ui^ denotes the advection 



U,j{x) = J (fk' M(fc')e-i'='-=«-5'^'(C'+c.)fc' (A5) 



and Sii denotes the deformation 



S.jix) = -iJ d2fc'M(fc')fc'^e-"^' ="-5'='(^-+"='^)'^' , (A6) 

which can be expressed as Uij{x)'^^ . We sort the ob- 
tained equations in orders of the V operator and obtain 



with Vi2 








Vi2 






/ 



Xi-V 1 Wu 

-af'-/3f^ \Wu. (A7) 

The matrix on the left side can be inverted and we obtain 
x'^-^J (A8) 

where we omitted terms proportional to V12 on the right- 
hand side. From eq. (A8) we obtain the evolution equa- 
tions 



N 



(A4b) 



Xi = ^TjUijixi - Xj) + q-^{xi - Xi{t)) , 

N (A9) 
Ci = H \Sij(xi — Xj)Ci + CiS^j(xi — 3;^)] 



(q + Q - 
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for the most probable evolution of elliptic vortices. 

An analogue calculation leads to the evolution equa- 
tions for the elliptic structures in the auxiliary field 



=0, 



N 



Xi — ^ ^ Fj IJij [xi Xj ) , 



(AlO) 



N 

Ci = ^ ^ Tj ^Sij {xi — Xj) Ci + CiSfj (xi -— Xj ^ 



where we defined the advection and deformation for the 
auxiliary field 

U,,{x) = [ d2fc'M(fc')e-"='— 

7 (All) 

S^,,{x)^U,,ix)^^. 



To derive these equations, one has to take into account 
the variation of the Lagrangian with respect to Ti, Xi 
and Ci as well as their time derivatives Fj, and Ci. 
Note, that we omitted higher order terms in Ci. 



Appendix B: Transition probability 



Sc = ^<^^kdtCj[—k,t)Q{k)(2j{k,t) where the k inte- 
gration is Gaussian and leads to 



S'o = Itt j At 



r? 



f2 
^ 2 



det(Q + 2Ci) ddei{Q + 2C2) 



Here we omitted the overlap in accordance with our as- 
sumptions. In the Instanton approximation Z w e''^" 
we get an zeroth-order approximation to the transition 
probability density 



Ic{{xin}\{x,,T,})^e 



iSe 



(Bl) 



which depends on the shape C of the vortices at 
time t* . The time t* has to be extracted from 
Ti{t) = qf,{t -t*)+Ti{t*). The transition probability 
/ can be written as 



f {{xlVl}\{x,,T,}) ^ 

jdCp{C)Jc{{xlT°,}\{x,,T,}) (B2) 



With the Instanton equation for the vorticity field 
u]{x,t) from equation Q we obtain the extremal action 



where p{C) is a distribution of initial shapes C . The 
investigation of this expression is left for future work. 
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